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Abstract 

The ground state and excited spectra of baryons containing three identical heavy quarks, b or c, 
have been recently calculated in nonperturbative lattice QCD. The energy of positive and negative 
parity excitations has been determined with high precision. Lattice results constitute a unique 
opportunity to learn about the quark-confinement mechanism as well as elucidating our knowledge 
about the nature of the strong force. We analyze the nonperturbative lattice QCD results by 
means of heavy-quark static potentials derived using SU(3) lattice QCD. We make use of different 
numerical techniques for the three-body problem. 

PACS numbers: 14.40.Lb,12.39.Pn,12.40.-y 
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I. INTRODUCTION 


A precise calculation of the mass of the ground state triply-bottom baryon, Qtbb, and ten 
excited positive and negative parity states with angular momentum up to J = 7/2 has been 
recently performed in nonperturbative lattice quantum chromodynamics (QCD) with 2 + 1 
flavors of dynamical light quarks in Refs. The b quark was implemented very accurately 


with improved lattice nonrelativistic QCD as suggested in Ref. |3[. There are also recent 
nonperturbative calculations of the ground and excited states with angular momentum up 
to J = 7/2 of triply-charm baryons using anisotropic lattice QCD with a background of 
2 + 1 dynamical light quark helds [^. Within the same framework of Refs. [H, |^, dynamical 
2 + 1 flavor lattice QCD and improved nonrelativistic lattice QCD for the b quarks, the 
bottomonium spectrum was computed in Ref. obtaining an excellent agreement with the 
experiment. A calculation of the charmonium spectrum with the same lattice action and the 
same lattice spacing used in Ref. can be found in Ref. [^, where one can get an idea of 
the typical size of the systematic uncertainties. Such studies grant the lattice QCD results 
for the bound state problem of three-heavy quarks with the category of experimental data 
and thus provide us with a unique opportunity to test phenomenological quark models for 
baryons in the energy regime where the description using potential models is expected to 
work best. 

As pointed out by Bjorken time ago [^, bound states of three-heavy quarks, QQQ, may 
provide a new window for the understanding of baryon structure. Baryons made of three- 
heavy quarks reveal a pure baryonic spectrum without light-quark complications and provide 
valuable insight into the quark-conhnement mechanism as well as to elucidate our knowledge 
about the nature of the strong force. On the theoretical side one would expect the potential 
models would be able to describe triply-heavy baryons to a similar level of precision as 
their success in charmonium and bottomonium. In the same way the QQ interactions are 
examined in heavy mesons, the study of triply-heavy baryons will probe the QQ interactions 
in the heavy quark sector. Triplwheavy baryons have been studied by different methods, 
including quark models (see Refs. M[^ for a comprehensive review of quark-model theoretical 


approaches), QCD sum rules [lCHl3l| and potential nonrelativistic QCD with static potentials 



from perturbation theory at leading order [1^ and next-to-next-to-leading order 
However, no experimental results are available so far for triply-heavy baryons (see Ref. 
for a recent calculation of production cross section at the LHC), and thus, the predictions 
of their properties cannot yet be compared to the real world. 

Our poor knowledge of the three-quark, 3Q, potential stems on the difficulty to produce 
QQQ states and, as said above, the consequent lack of experimental data. The precise 
calculation of the ground and excited states of triply-bottom baryons 0,0, together with 
the ground and excited states of triply-charm baryons 0, provide us with the ideal testbench 
to improve our understanding of phenomenological quark models for baryons in the heavy- 
quark sector. The quark-model dependent calculations could be tested by comparing them 
to nonperturbative hrst-principles calculations in lattice QCD of the bbb and ccc systems. An 
adequate starting point may be the static 3Q potential derived in SU(3) lattice QCD that 
will be described in the next section. In Sect. lIHI we will present our results making emphasis 
in the comparison between nonperturbative lattice QCD data and the results obtained with 
the static 3Q potential derived in SU(3) lattice QCD. For this purpose we will make use of 
different numerical techniques for the three-body problem, generalized Gaussians variational 


approaches, hyperspherical harmonics and Faddeev equations |18l-l20[|. to test quark-quark 
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potential descriptions in the heavy quark limit. Finally in Sect. IIVI we will briefly summarize 
the main conclusions of this study. 


II. THE STATIC THREE-QUARK POTENTIAL IN SU(3) LATTICE QCD 


Since the early days of QCD the interaction among heavy quarks has been explored as an 
important tool to learn about the behavior of QCD at low-energies. The QQ static potential 
has been studied extensively by lattice gauge theories [2l[, being nowadays a very well-known 
quantity. The typical shape of the color-singlet QQ static potential, which is characterized 
by a short-range Coulomb behavior and a long-range linear rise, well represents the double 
nature of QCD as an asymptotically free and infrared confined theory. The excitation 
spectrum of the gluon held around a static quark-antiquark pair has also been explored by 
lattice calculations jl^. On the large length scale the spectrum agrees with that expected 
for stringlike excitations while in the short range it shows a Coulomb-like behavior as it was 
first noted within the context of the static bag picture of gluon excitations j^ . 

In QCD the three-quark potential is of prime importance reflecting the SU(3) gauge 
symmetry and being directly responsible for the structure and properties of baryons, in 
the same way the QQ potential is responsible for the meson properties. Furthermore, the 
3Q potential is a key quantity to clarify the quark confinement mechanism in baryons. 
However, up to now lattice QCD studies have paid relatively less attention to the potential 
that describes the interaction of three-heavy quarks, as a consequence of their still missing 
experimental evidence. Thus, the 3Q potential is much less known than the heavy QQ 
potential for which many lattice studies exist j2l|. Even for the ground-state, its accurate 
measurement using lattice QCD is relatively difficult and has been performed recently 


27| . The gluonic excitation (excitations of the gluon field that are not thus of quark origin) 

concluding 


of the three-quark system has also been studied in SU(3) lattice QCD [2^, [29 
that the lowest excitation would be almost 1 GeV above the ground state. This is rather 
large in comparison with the low-lying excitation energy of the quark origin what makes 
the gluonic excitation mode invisible in the low-lying excitations of baryons. Thus, quark- 
degrees of freedom play the dominant role in low-lying baryons with excitation energy below 
1 GeV. This large gluonic excitation energy is conjectured to give a physical reason of the 
success of the quark model for low-lying baryons even without explicit gluonic modes j^ . 

Most of the existing lattice studies of the three-quark static potential have mainly ex¬ 
plored the region of large interquark distances j24-3^. As for the QQ case, the characteristic 
signature of the long-range non-Abelian dynamics is believed to be a linear rising of the static 
interaction. Moreover, the general expectation for the baryonic case is that, at least clas¬ 
sically, the strings meet at the so-called Fermat (or Torricelli) point, which has minimum 
distance from the three sources (V—shape configuration). If this is the case, one should see 
a genuine three-body interaction among the static quarks. The confining short-range 3Q 
potential could be also described as the sum of two-body potentials (A—shape configura¬ 
tion) 0, im, il' 

Y— from the A- 


M. Many of the lattice calculations have focused on distinguishing the 
configuration, despite the difference between a A and a V—shape potential 
being rather small and difficult to detect. 

In recent years, according to the remarkable progress of the computational power, the 
lattice QGD Monte Garlo calculations have become a reliable and useful method for the 
analysis of nonperturbative QGD. In particular, the QQ potential, responsible for the meson 
properties, has been extensively studied using lattice QGD. The data of the QQ ground state 
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TABLE I: Standard string tension, a, Coulomb coefficient, A, and constant term, C, of the static 
heavy quark potential obtained from SU(3) lattice QCD taken from Ref. 2^ for a lattice spacing 
a = 0.19 fm. 



cj (a 2) 

(T (GeVQ 

A 

C (a-Q 

C (GeV) 

QQ 

0.1629(47) 

0.1757 

0.2793(116) 

0.6203(161) 

0.6442 

3Q(Y) 

0.1524(28) 

0.1644 

0.1331(66) 

0.9182(213) 

0.9536 

3Q(A) 

0.0858(16) 

0.0925 

0.1410(64) 

0.9334(210) 

0.9694 


potential are well reproduced by a sum of a Coulomb term, due to the perturbative one-gluon 


exchange (OGE), and a linear conhnement |25l. l29 


yQQ(r) = 


^QQ 


+ <^(3(3 + Coo 


QQ 


( 1 ) 


We give in Table |I]the parameters of the QQ potential taken from Ref. 2^. The value of the 


QQ conhnement strength was determined to reproduce the value obtained from the linear 
Regge trajectories of the pseudoscalar tt and K mesons, ^/a = (429±2) MeV (^ . 

However, there is almost no reliable formula to describe the 3Q potential directly based 
on QCD, in spite of its importance for the study of baryon properties. There have been 
recent advances on the determination of the ground state 3Q potential that is expected to 
take the form 251, 


) - -^3Q X] Tvr^- 


rA' T 
*^3(3 -^min 


+ C[ 


3(3 5 


( 2 ) 


i<j 


where Tmin is the minimal value of the total length of color hux tubes linking the three 
quarks, and 


^ 3(3 — —Ann 


QQ 1 


^Iq — ^QQ ■ 


(3) 


The short-distance behavior of V^Q^r) is expected to be described by the two-body Coulomb 
potential as the one-gluon exchange result in perturbative QCD. Fit analysis in terms of 
a Yukawa potential considering a possible gluon mass have concluded = 0, what 
reduces the Yukawa to the Coulomb potential jl^. The one-gluon exchange result indicates 
also the simple relation on the Coulomb coefficients in the QQ and the QQQ potentials as 

A 3 Q z 

like the Y—ansatz 29, ^ 1. 


~ ^Aqq. For tlm long-distance behavior, the confining baryonic static potential rises 


Calculations both in full and quenched QCD demonstrates that the conhning baryonic 
static potential approaches the A—ansatz at short distances 0 H ij- The A —ansatz 
behavior at short distances is of great importance for phenomenological models since the 
calculation of orbital excited states with the A—ansatz are much more simpler. For an 
equilateral QQQ arrangement, as expected for a system of three-heaw identical quarks, 
departure from the A—ansatz is not signihcant until dgg ~ 0.7 fm 1271 . l32l . so that the 


A—ansatz may be the more relevant one for quarks conhned inside a baryon and then, for 
the case of QQQ bound states whose root-mean square radii are much smaller than such 
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distance (see Table HB . Thus, the static three-quark potential could be described by a simple 
sum of the effective two-body QQ potentials with a reduced string tension j25l-l27l| , 




3Q 


r = 


-^3Q ^ ^ 


a. 


3Q 


i<j 


where 


A3Q ~ - A 


QQ y 


(To 



(4) 

i<j 


0.53(Tgg . 

(5) 


The reduction factor in the string tension can be naturally understood as a geometrical 
factor rather than the color factor, due to the ratio between and the perimeter length 
of the 3Q triangle, suggesting = (0.50 ~ 0.58) (T qq [ 2 ^. For the particular case of quarks 
in an equilateral triangle = 0.58(Tgg |32[. When the A— ansatz is adopted for 

the two-body linear potential, still the same relation holds for the strength of the Coulomb 
potential A^q ~ |Agg. This A— ansatz has been widely adopted in the nonrelativistic quark 
model because of its simplicity SB [3 [m, 


III. RESULTS AND DISCUSSION 


In order to check if the static three-quark potential of Eq. (|1]) (or Eq. ([2])) with the 
parameters in Table [I] determined from lattice QCD can reproduce the bbb and ccc baryon 
spectra measured also in lattice QCD S, B 0] we will make use of three different numerical 
methods: generalized Gaussians variational approach jl^, hyperspherical harmonics [l9[ 
and Faddeev equations [l8|. The three-methods have been used and the difference in results 
is negligible. In all cases we solve the nonrelativistic Schrodinger equation 


{Ho + V,Q{r)}^ir}=E'^ir), 

where Ho is the free part of three-heavy quarks without center of mass motion 


/fo = ^(Mg + 


i=l 


Pi 


2M< 


Q 


T, 


CM 


and Mq is the mass of the heavy quark. The mass of the heavy baryon will be finally given 
by = 3Mq -|- E. 

We have checked the validity of the A—ansatz for the spectroscopy of triply-heavy baryons 
calculating their masses and root mean square radii by means of a simple confining interac¬ 
tion given either by a U—shape or a A—shape interaction. The U—shape potential would 
be given by j^ . 

^3q(t) = CT^g Emin , (6) 

where 


Lmin 


1 /o 

- (a^ + b'^ + c^) -h —\/{a + b + c){-a + b + c){a-b + c){a + b-c) 


1/2 


(7) 


being a, b and c the three sides of the 3Q triangle, when the angles of the 3Q triangle do 
not exceed 2tt/3. When an angle of the 3Q triangle exceeds 27r/3, one gets, 

Lmin = a + b + c — max (a, b, c) . (8) 
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TABLE II: Mass, in GeV, and root-mean square radius, in fm, of the ground state, and 

the first positive parity excited state, £’2(1'*'), of QQQ baryons for different values of the mass of 
the heavy quark Mq, in GeV, and for the V—shape and A—shape confining potentials of Eqs. ([6]) 
and ([9]), respectively. 


Mq 

Y- 

^i(i^) 

-shape 

^2(r) 

^i(i^) 

A- 

V<ri> 

-shape 

A 

V 

> 

4.0 

12.814 

13.121 

12.795 

0.236 

13.094 

0.334 

4.4 

13.989 

14.286 

13.970 

0.228 

14.260 

0.324 

4.8 

15.166 

15.455 

15.148 

0.222 

15.430 

0.315 

5.2 

16.346 

16.627 

16.328 

0.216 

16.603 

0.306 

5.6 

17.528 

17.801 

17.509 

0.211 

17.778 

0.299 


We have taken for the string tension a standard valne of o^q = 0.1648 GeV^ 

A—shape potential wonld have the form, 

IAq( r) = X] |g - fj I, (9) 

i<j 

with a^Q = 0.53(t|^ = 0.0874 GeV^. The resnlts are shown in Table UTl for different valnes 
of the heavy qnark mass. As can be seen the two geometrical confignrations give the same 
resnlt with a small difference of 0.15%. If the geometrical factor in front of the A—shape 
interaction is slightly modihed, 0.55 instead of 0.53, the V—shape and the A—shape resnlts 
wonld be exactly the same. This result was already noted in Ref. ji^ . Besides, one observes 
how the root-mean square radii are much smaller than the quark-quark distance for which 
the Y— and A—shapes start to slightly differ, ~ 0.7 fm j^ . 

Looking for the general pattern of the results of nonperturbative lattice QGD for the bbb 
system we have compared against the result of a single A—shape conhning potential or a 
single Goulomb interaction. The results are shown in Fig. [1] In both cases we have hxed the 
strength of the potential to reproduce the £ 2 ( 3 / 2 +) — £i(3/2+) mass difference correctly, 
obtaining = 0.2076 GeV^ for the case of the single A—shape conhning interaction and 
A^q = 0.4231 for the single Goulomb potential case. We clearly note the large strength 
of the two interactions as compared to the predictions of SU(3) lattice QGD (see Tabled]), 
where both terms would simultaneously contribute. The A—shape interaction (upper panel) 
produces a splitting between the positive (i.e., the hrst excited state of = 3/2+) and 
negative (i.e., the hrst excited state of = 3/2“) parity excited states too large and 
a rather small excitation energy for the positive parity excitations. On the contrary, the 
Goulomb interaction (lower panel), being an almost hypercentral potential, drives to the 
expected degeneracy between the positive and negative parity excited states, not observed 
in the nonperturbative lattice QGD results. Thus, as the strength of the potential has been 
hxed to reproduce the £ 2 ( 3 / 2 +) — £i(3/2+) mass diherence correctly, the negative parity 
excited states are obtained close to the positive parity excited states in clear disagreement 
with the nonperturbative lattice QGD results. Besides, the positive parity states (except 
the hrst excited state of = 3/2+ that has been htted) are predicted slightly above the 
nonperturbative lattice QGD results. 


24 -29 . The 
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FIG. 1: bbb excited state spectra, solid lines, for a single A—shape confining potential (upper 
panel) or a single Coulomb interaction (lower panel). In both cases we have fixed the strength of 
the potential to reproduce the E2{3/2'^) — Ei{3/2'^) mass difference correctly, obtaining a^g = 
0.2076 GeV^ for the case of the single A—shape confining interaction and A^q = 0.4231 for the 
case of the single Coulomb potential. The boxes stand for the nonperturbative lattice QCD results 
of Ref. 13. 


Thus, bbb nonperturbative lattice QCD calculations point to a static 3Q potential given 
by a mixture of a A—shape confinement and a Coulomb interaction. We have therefore 
make use of the standard 3Q static potential derived from SU(3) lattice QCD and shown in 
Eq. (jl]) with the parameters reported in Ref. [l^] and quoted in Table [H A^q = 0.1410 and 
a^Q = 0.0925 GeV^ to calculate the bbb and ccc spectra. The results are shown in Fig. [2l We 
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FIG. 2: 656 (upper panel) and ccc (lower panel) excited state spectra, solid lines, for the 3Q static 
SU(3) lattice QCD potential of Eq. (j4]) with the parameters of Tabled! The boxes stand for the 
nonperturbative lattice QCD results of Ref. Q] for the 666 system and Refs, 
system. 


4ll | for the ccc 


have fixed the mass of the heavy ^ark to reproduce the = 3/2+ ground state of the bbb 
and ccc systems, 14.372 GeV [ll, lH and 4.758 GeV [^, respectively, obtaining mb = 4.655 
GeV and me = 1.269 GeV. As we can see in Fig. [2l there is a large difference in the excited 
states with the nonperturbative lattice QGD results, predicting a small splitting between 
positive and negative parity excited states and also a small excitation energy for the positive 
parity states. These results clearly point to a lack of strength either in the confining or in 
the Goulomb potential. 
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One may think about several possible reasons for the disagreement: 

• The parameters of Eq. dH) obtained in quenched QCD in Ref. [1^ might be different 
from those in 2+1 flavor QCD employed in Refs. 0,31. To check this possibility, 
one has to calculate parameters of Eq. (jl]) in 2+1 flavor QCD but such calculations 
do not exist so far. 

• Results from both Ref. and Refs. 0,33 contain systematic errors such as lattice 
artifact and chiral extrapolation as well as statistical errors. This may cause the 
disagreement. 

• The fitting form Eq. (jTj), or Eq. ([2]), might not be appropriate to describe the static 
three-quark potential in lattice QCD. 

• The static quark description might be inaccurate for bbb and, in particular, for ccc. 
Higher order terms might be necessary. 

• The quark model description with ” 3-quark potential” might be inappropriate for bbb 
and ccc systems. 


IV. SUMMARY 


To summarize, the spectra of baryons containing three identical heavy quarks, b or c, 
have been recently calculated in nonperturbative lattice QCD. The energy of the lowest 
positive and negative parity excited states has been determined with high precision. These 
achievements constitute a unique opportunity to test phenomenological potential models in 
the regime where they are expected to work best. We have analyzed these results by means 
of static three-quark potentials derived using SU(3) lattice QCD using different numerical 
techniques for the three-body problem. Our results confirm the expectations of SU(3) lattice 
QCD of an almost indistinguishable confining Y— or A—type at short distances for heavy 
baryons. The static three-quark potential with parameters from lattice QCD does not 
reproduce bbb and ccc spectra. 

At the light of the present results a further effort to obtain a constituent quark-model 
potential description of the nonperturbative lattice QCD results, as has been done in the 
past for the heavy-meson systems 


-3, may help in understanding the connection between 
static three-quark potential parameters and simple Cornell-like potential descriptions. This 
work is in progress 
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